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We present a comprehensive experimental and theoretical investigation of the thermodynamic 
properties: specific heat, magnetization and thermal expansion in the vicinity of the field-induced 
quantum critical point (QCP) around the lower critical field H cl « 2 T in NiCl 2 -4SC(NH 2 ) 2 . A T 3/2 
behavior in the specific heat and magnetization is observed at very low temperatures at H = H c i that 
is consistent with the universality class of Bose-Einstein condensation of magnons. The temperature 
dependence of the thermal expansion coefficient at H c \ shows minor deviations from the expected 
T 1//2 behavior. Our experimental study is complemented by analytical calculations and Quantum 
Monte Carlo simulations, which reproduce nicely the measured quantities. We analyze the thermal 
and the magnetic Griineisen parameters that are ideal quantities to identify QCPs. Both parameters 
diverge at H c i with the expected T~ power law. By using the Ehrenfest relations at the second order 
phase transition, we are able to estimate the pressure dependencies of the characteristic temperature 
and field scales. 

PACS numbers: 75.40.-s; 64.70.Tg; 65.40.De 



I. INTRODUCTION 

Bose-Einstein condensation (BEC) has triggered great 
interest in the last years and was found in a variety of 
complex many body systems, such as cold atoms, su- 
pcrfluid Helium or superconductors. By exploiting the 
Matusbara-Matsuda mapping of S = 1/2 spins into hard 
core bosons PBatyev showed that the field induced phase 
transition between canted XY antiferromagnetic (AFM) 
ordering and the fully polarized state can also be de- 
scribed as a BECP This useful mapping between mag- 
netic systems and dilute gases of bosons can be extended 
to higher spin values^' and it has been successfully ex- 
ploited first on TlCuCl^ and other quantum magnet sP^ 

One material investigated recently, is MCI2- 
4SC(NHa)2P also known as dichlorotetrakisthiourea- 
nickel (DTN). It has a body-centered tetragonal crystal 
structure with chains of Ni-Cl-Cl-Ni atoms arranged 
along the crystallographic c direction. DTN enters the 
XY-AFM ordered state between moderate fields of 2T 
and 12.5 T, if the magnetic field H is applied along c. 
The magnetic atom in DTN is Ni 2+ carrying a spin 
S = 1 due to an almost completely quenched orbital 
momentum. The Hamiltonian for DTN can be written 
as 

H = J2 JvS, • S r+e „ + ^[£(^) 2 - guBHS*}, (1) 

ru r 

where v = {a, b, c} and e„ is the relative vector between 
nearest-neighbors along the ^-direction. The magnitude 



of the dominant single-ion anisotropy D is 8.9 KP The 
AFM exchange interactions between neighboring spins 
are J c — 2.2 K along the chains, and about 10 times 
smaller, J a b =0.18K, in the a6-plane. The last Zeeman 
term in equation |l]) is originated from the applied mag- 
netic field H and the quantization z-axis is chosen along 
the field direction. The gyromagnetic factor g parallel to 
the c-axis was estimated to be 2.26 by ESR experimentsP 
Equation ([!]) is only a minimal Hamiltonian for describ- 
ing the magnetic properties of DTN. Further contribu- 
tions, such as dipolar interactions, which break the U(l) 
symmetry of global rotations along the spin z-axis are 
small, but they become relevant at very low tempera- 
tures. Therefore, the critical exponents characteristic of 
a BEC quantum critical point (QCP) can only be ob- 
served, if the U(l) symmetry breaking terms are at least 
one order of magnitude smaller than J a t. At low enough 
temperatures one should observe a crossover from the be- 
havior characteristic for BEC-QCPs to the one expected 
for an Ising-likc QCP. Investigations of the exact shape 
of the phase boundary close to H c i and H c2 down to 
1 mK by detailed AC susceptibility measurements evi- 
denced the universality class of a BEC in DTNP^ Up to 
date, this is the solely experimental observation consis- 
tent with a field induced BEC-QCP in this material. 

The universality class of the QCP can also be deter- 
mined by measuring the exponents for the power-law de- 
pendencies of different thermodynamic quantities as a 
function of temperature. Table|l] shows the expected ex- 
ponents for BEC and Ising-like QCPs in 2 and 3 dimen- 
sions dP It is important to note, that d = 3 is the upper 
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TABLE I: Temperature dependencies 
namic quantities: magnetization M(T), 



of the thermody- 
thermal expansion 



a(T) and specific heat C(T) at the field induced QCP. 
The variable d denotes the spatial dimensionality of the sys- 
tem. The exponents of the Ising-like QCP are given for d — 3. 



critical dimension for the Ising-like QCP (D = d+z = 4), 
where z = 1 is the dynamical exponent and D the effec- 
tive dimensionality. Therefore, one should expect further 
logarithmic corrections to the power law behaviors listed 
in Table]!] In this work, an extensive study of the mag- 
netization M(H,T), specific heat C(H,T), and thermal 
expansion av(H,T) close to the critical field H c \ gives 
further strong evidence that DTN belongs to the univer- 
sality class of BEC. 

It has been shown recently,^ that QCPs can be 
detected by measuring the divergence of the thermal 
Griineisen parameter, 



ay 
~C 



(2) 



for pressure tuning, and the magnetic Griineisen param- 
eter, 



dM/dT 



(3) 



for a magnetic field tuned QCP. Both Griineisen param- 
eters diverge at the QCP like T oc T~ x l vz , where v is the 
critical exponent that relates the correlation length with 
the driving parameter of the quantum phase transition 
(magnetic field in the case of DTN). The effective dimen- 
sion D — d + z is higher or equal to 4 with z = 1 for the 
Ising-type and z = 2 for the BEC-type QCP, and we get 
v = 1/2 for both cases. Therefore, Y oc T _1 is expected 
for a BEC-QCP and Y oc T~ 2 for an Ising-like QCP in 
agreement with the power laws listed in Table]!] 

Typical for QCPs is the occurrence of local maxima 
in the entropy due to enhanced quantum critical fluctu- 
ations. This implies a sign change of the thermal ex- 
pansion coefficient which is linked to the entropy via 
the Maxwell relation ay — —dS/dp. In the past, the 
concept of the diverging Griineisen parameter was used 
successfully to identify and characterize not only well un- 
derstood magnetic QCPs, but also other more puzzling 
QCPs in intermetallic compounds.^ 

So far, dilatometric properties were used to investi- 
gate the quantum critical behavior of only a few quan- 
tum magnets. The d = 3 coupled spin-dimer system 
TlCuCLP and the quasi-one-dimensional spin- ladder 



compound (C5Hi2N)2CuBr4i^ are two rare examples. 
These system s show field induced phase transitions at low 
temperatures pESEZI however, investigations of the ther- 
mal expansion coefficient a and the Griineisen ratio Y t h 
show significant deviations from the expected behavior 
in both cases. Dilatometric experiments on TICUCI3 re- 
veal that while the thermal Griineisen parameter diverges 
with the expected power law 1/T, the individual quanti- 
ties specific heat C and thermal expansion a fail to fol- 
low the predictions.^ In (CsH^N^CuBr-i, the thermal 
expansion coefficient a c along the crystallographic c di- 
rection shows a weak indication of 1/yT divergency and 
a clear sign change at the lower and upper critical field 
H c i/c2, but a detailed discussion of Y t h is missing.^ The 
present study of the thermal and magnetic Griineisen ra- 
tios shows that DTN is an excellent candidate to close 
this gap of knowledge about dilatometric properties of 
insulating quantum critical materials. 

The paper is organized as follows: in Section]!!] we de- 
scribe the experimental techniques that we used in static 
magnetic fields to measure the specific heat, thermal ex- 
pansion and magnetization up to 15 T, the approximated 
model used for the analytical calculations, and the quan- 
tum Monte Carlo (QMC) simulations of the thermody- 
namic quantities. Section [TTT] contains a detailed descrip- 
tion of the experimental and theoretical results. We con- 
tinue in Section [TV] with the comparison between experi- 
ment and theory and the analysis of the anomalies at the 
phase boundary via the Ehrenfest relations. Section [V] 
summarizes the most important results of our study. 



II. METHODS 

The preparation of high-quality single-crystals is ex- 
plained elsewherePSl All experiments were conducted, 
partly down to 30 mK, inside commercial available dilu- 
tion refrigerators, furbished with superconducting (SC) 
magnets with maximum fields of up to 15 T. We mea- 
sured the magnetization with a high resolution Fara- 
day magnetometer.^] The thermal expansion and mag- 
netostriction experiments were carried out with a high 
precision capacitive dilatometer^l made of CuBe. The 
dilatometer can be rotated by 90° in order to measure 
the length change not only parallel but also perpendicu- 
lar to the applied magnetic field. The specific heat was 
measured with the compensated heat pulse technique^ 
and the data were confirmed by experiments using the 
dual-slope method^ on the same sample platform. The 
precise match between both sets of experimental data is 
remarkable. Additionally, we used the specific-heat setup 
to perform magnetocaloric effect (MCE) measurements 
for a precise estimate of the critical field H c \ . 

The analytical calculations of the various thermody- 
namic properties were based on the usual expansion in 
the gas parameter or ratio between the scattering am- 
plitude and the average inter-particle distance, p -1 / 3 |2H 
For this purpose we mapped the S z = 1 magnetic excita- 
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tions of the low field paramagnetic (PM) state into hard 
core bosons, where the z-component of the magnetiza- 
tion density in the original model, (S*), is mapped into 
the particle density p. Here we neglect the contribution 
to the magnetization of the S z — — 1 modes, because we 
are assuming that H is close to H c \ and T -C A, with 
A ~ 3K being the H = spin gap of DTN. We use 
the expression derived in Ref. [S7J for the single-particle 



dispersion , w k 



gf-isH with 



s 2 e k . 



(4) 



The parameters s 2 and /i are given by the following ex- 
pressions: 



s z = 2- 



1 

N 



E 



fx + s-e k 



k K 



By using the Hamiltonian parameters for DTN esti- 
mated in Ref. [S] the resulting values are s 2 = 0.92 and 
/i = 10.3 K. The effective repulsion between bosons in the 
long wavelength limit, vq = r (Q,Q) with Q being the 
ordering wave- vector, results from summing the ladder 
diagrams for the bare interaction vertex Vy 2 ^ 



r q (k,k') = v q - r 

J —■ 



d P 3 r p (k,k / ) 

87T 3 Wk+p + W k /_ p 



(6) 



where V q = C/ + 2 J c cos q z + 2 J a (, (cos g K + cos q y ) for DTN 
and U — > oo is included to enforce the hard core con- 
straint. The effective Hamiltonian in the long wavelength 
limit |k — Q| <C 1 is given by 



H 



ell' 



k 



(ek-Ai)a k a k + 



2N 



E 

k.k'.q 



«k+q a L 



-q a k a k'' 



(7) 



where TV is the total number of lattice sites and the opera- 
tor a k (a k ) creates (annihilates) a boson with momentum 
k. e k is obtained by taking the long wavelength limit of 
w k : 



2m* 



[k x + k y ) 
2m* 



(8) 



The resulting quadratic mean-field Hamiltonian can be 
easily diagonalized and the various thermodynamic prop- 
erties are computed by solving the self-consistent condi- 
tion imposed by Eqs.(10l and (11) 

The analytic calculations have been supplemented 
by large scale numerical simulations of the microscopic 
model. We have used the Stochastic Series Expansion 
(SSE) QMC method to simulate the Hamiltonian (1) on 
finite-sized lattices using the experimentally determined 
parameters. The SSE is a finite-temperature QMC tech- 
nique based on importance sampling of the diagonal ma- 
trix elements of the density matrix e~@ H POUSI The use 
of operator loop cluster updates reduces the autocorrela- 
tion time for the system sizes. We consider here up to 
«2x 10 4 spins to at most a few Monte Carlo sweeps even 
at the critical temperature.^ This enables us to explore 
the vicinity of the critical points very efficiently. On the 
dense temperature grids needed to study the critical re- 
gion in detail, the statistics of the Monte Carlo results 
can be significantly improved by the use of a parallel 
tempering scheme!- 31 * 33 ! The implementation of this tem- 
pering scheme in the context of the SSE method has been 
discussed in detail previously) 32 * 34 ! and we follow the one 
developed in Ref. [51]. 

QMC estimates for observables of a spatially 
anisotropic system can depend non-monotonically on the 
system size for isotropic lattices. One can instead use 
anisotropic lattices to more rapidly obtain monotonic 
behavior of the numerical results for extrapolating to 
the thermodynamic limit. Anticipating similar effects in 
the present model (since J a , Jb *C J c ), we have studied 
tetragonal lattices with L x = L y = L z /4. 

The specific heat has been extracted from the simu- 
lation data by the numerical differentiation of the total 
internal energy of the system - a quantity that is esti- 
mated extremely accurately by the SSE method. The 
temperature dependence of the energy is approximated 
by a polynomial in T, and the derivative of the poly- 
nomial fit is used to estimate the specific heat. Thus, 
artifacts of discrete numerical differentiation of the raw 
data are avoided and we yield a relatively noise-free spe- 
cific heat curve. 



with 



<9 2 cj k 



dkl 



(9) 



k=Q 



The chemical potential fi is gpsH — wq. After a mean- 
field treatment of W c g in the PM phase, H < H c i , the 
interaction term simply leads to a renormalization of the 
chemical potential, fi — > p,, with 



A = M - 2v oP 
and the particle density 

P= ^E< a k a k>- 



(10) 



(11) 



III. EXPERIMENTAL RESULTS 

A. Magnetization 

Figure[T] shows a comparison of the magnetization M 
as a function of magnetic field measured perpendicular 
at T = 500 mKP and parallel at T = 16 mK^ to the 
crystallographic c-axis. The latter data nicely reflect the 
Neel-ordered state in the ab plane with increasing canting 
along c between 2.1 T and 12.6 T followed by saturation. 
For H || ab, the magnetization shows PM behavior with 
no ordering, approximately following a Brillouin function 
and saturating at 2.2 hb per Ni 2+ -atom around 15 T. In 
this field direction the magnetic field increases the size of 
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FIG. 1: (color online) Magnetization M as a function of the 
magnetic field H of DTN at 16 mK for measurements H par- 
allel (solid line, Ref. 20 ) and perpendicular to the crystal- 
lographic c direction at 500 mK (dashed line, ref. [H]). The 
inset shows new results of the magnetization at 50 mK in low 
fields up to 2.5 T and a linear fit (dashed line) to the data. 
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FIG. 2: (color online) Corrected magnetization AM versus 
T 3/2 for magnetic fields H || c between IT and 3T. The 
magnetization at the critical field fJ,oH C 2 = 2.08 T is labeled 
with stars. Data above IT are shifted vertically by 0.02 /is 
per Ni atom for better visualization. Arrows indicate the 
phase transition into the XY-AFM state. 



the spin gap instead of closing it as happens for H || c. 
The inset of Fig.[T] shows the low field part of the mag- 
netization H || c at 50 mK in greater detail. We observe 
a linear increase between zero and the lower critical field 
H c i, which can not be explained with a U(l) invariant 
Hamiltonian such as T~L, where the magnetization is sup- 
posed to be zero at T = in the quantum PM region 
H < H c i. This effect cannot be caused by single ion 



impurities, e.g. uncoupled Ni 2+ moments, because such 
spins should be fully polarized for magnetic fields well 
below the lower critical field H -C H cl . Instead, we con- 
clude that the linear slope (dotted line) is caused by a 
misalignment of the sample, which gives a contribution 
M(H _L c) to the magnetization. From the value of the 
susceptibility compared to data for H _L c we estimate a 
misalignment of less than 1.2°. 

Furthermore, we measured the magnetization for H \\ c 
near H c \ and in the temperature range 0.1 K< T <0.5 K 
(data not shown) to extrapolate the phase boundary for 
T -> 0. 10 We obtain a critical field H cl of 2.08 T. It is 
important to note that this value depends on the specific 
conditions under which the sample is mounted in the ex- 
perimental setup. The slightly higher value of H c \ in the 
magnetization compared to the specific heat and thermal 
expansion L || c values (see data below) supports the as- 
sumption of small sample misalignment, because angu- 
lar resolved measurements of the magnetostriction have 
shown that H c \ increases with increasing angle between 
the field direction and the crystallographic c-axisf^ 

The small PM contribution Mpm to the magnetization 
made it necessary to refine the M(T) data as a function 
of temperature for H \\ c. Figure[2] shows the corrected 
values AM(T) = M - M PM plotted versus T 3 / 2 for a 
variety of different magnetic fields 1 T < [i^H <3T. The 
magnetization is exponentially suppressed at low temper- 
atures inside the quantum PM state (noH = 1 T) and 
develops a T 3 / 2 -behavior when approaching the critical 
field H c i = 2.08 T in agreement with the expected behav- 
ior for a BEC-QCP in 3 dimensions. The entrance into 
the XY-AFM state well above H cl is marked with a dip 
in the magnetization (arrows in Fig. [2]). Below the mini- 
mum inside the AFM phase, the magnetization increases 
with a power law for decreasing temperature. 



B. Specific Heat 

We estimate the critical field H c i by MCE measure- 
ments (data not shown) following the analysis of the 
highest slope (dT/dH) of the temperature during field 
scans^ and find a value H c \ — 2.06 T in the specific-heat 
experimental setup. 

The specific heat curve of DTN contains three contri- 
butions: nuclear Schottky (NS), magnetic Schottky (MS) 
and quantum-critical (QC) contributions. Each of them 
dominates in different regions of the H — T phase dia- 
gram. The specific heat caused by phonons can be ne- 
glected in the temperature range below 5 K. The NS con- 
tribution is difficult to master in this material, because 
it originates from several nuclei (H, N, CI) generating 
a huge fraction of the specific heat at temperatures be- 
low 0.1K, see e.g. in Fig. [3] the measurement at 2.06 T. 
In addition, the effective magnetic field (sum of the ex- 
ternal field and the field generated by the ordered mo- 
ments) becomes rather high in the AFM ordered state 
above H c \- This increases the splitting of the nuclear 



5 




1 0.1 
7(K) 7-(K) 



FIG. 3: (color online) Main panel: total specific heat C/T 
as a function of temperature T at the critical field H c \ = 
2.06 T (filled dots) and the resulting Cmag/T (open squares) 
after subtraction of the nuclear Schottky (NS) contribution 
at low temperatures. Inset (a) shows the data of the main 
panel as C/T versus T -3 to illustrate the NS specific heat as 
indicated by a dashed line. The prefactor of the NS specific 
heat Cns/T = aT~ 3 is plotted in inset (b) versus H 2 . It does 
not follow the expected H 2 field dependence. 



FIG. 4: (color online) Magnetic specific heat C ma g/T as a 
function of temperature T for magnetic fields between and 
2.2 T in a double-log scale. The arrows indicate the AFM 
phase transition for fields H > H c \. The critical field ob- 
tained from MCE experiments is 2.06 T. The inset shows the 
energy gap between and 1.7T, estimated from the exponen- 
tial temperature increase of the magnetic specific heat at low 
temperatures. 



energy levels further and the Schottky anomaly becomes 
very big. The inset (a) of Fig. [^demonstrates the way we 
subtracted the NS contribution from the original specific 
heat at the critical field /io-ff = 2.06 T. We fit the data 
between 80 mK and 50 mK as C/T = aT~ 3 , which is a 
good approximation for the high temperature behavior 
of the Schottky anomalyP^ The prefactor a for the NS 
specific heat is given as a function of magnetic field in 
the inset (b) of Fig.[3j The pre-factors for different nu- 
clei should be additive, a = J2 a i> f° r the case that all 
the different nuclear energy levels are in the high tem- 
perature limit (ksT 3> A„„ c / ear ). Because the Zeeman 
splitting is linear in magnetic field, a(H) should obey a 
H 2 dependence for H < H c \. This is not observed in 
DTN and presently not understood. 

The insulating behavior of DTN is challenging for spe- 
cific heat experiments, because no free electrons con- 
tribute to the thermal conductivity. Heat is carried only 
by magnetic excitations and by phonons, whereas the 
phonon contribution is negligibly small in the temper- 
ature range below 1 K. Therefore the different thermo- 
dynamic subsystems (nuclear spins, magnetic moments, 
lattice) are only poorly coupled to each other at low tem- 
peratures. This causes an out-of-equilibrium state of the 
sample at very low temperatures. Thus, we only include 
data above 80 mK in our analysis. 

Figure]!] shows the magnetic specific heat C mag /T = 
C/T — Cns/T in a double-logarithmic display after the 
subtraction of the nuclear Schottky contribution, Cjvs, 
between and 2.2 T. The broad maximum around 2.5 K 



in the zero field measurement is caused by the thermal 
population of the \S Z = ±1 > excited states that are a 
rather broad band due to dispersion caused by the ex- 
change interactions.^ This dispersion is also responsible 
for the only slightly shift of the maximum to lower tem- 
peratures with higher magnetic fields. For small fields, 
the specific heat data can be fitted with an exponential 
function, exp(— ^Ap) in the low temperature limit allow- 
ing us to extract the spin gap A. The inset of Fig. [4] 
shows the gap values, estimated from the experimental 
data. They decrease linearly from 3.22 K down to 1 K 
when H varies between 0T and 1.7 T, whereas the zero- 
field value is in close accordance with previous suscep- 
tibility results of 3.3 Kpl Below 1.7 T, the temperature 
range for exponential behavior is too small for reliable 
data fitting. From the zero field gap, A = 3.22 K, we can 
estimate the critical field, 



ti c \ — , 

9Vb 



where the gap closes (kg — 1.380 • 10 23 JK 1 and hb — 
9.274 • lO^JT- 1 ). The calculated value H cl =2.12T 
matches the experimental values summarized in Table[TT| 
within 5%. The arrows in Fig. [J] indicate the anomalies 
caused by the phase transition into the XY-AFM state 
in the 2.14 T and 2.2 T measurement. 
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FIG. 5: (color online) Linear thermal expansion coefficient a c 
measured along the crystallographic c-axis (H || c || AL C ) as 
a function of the temperature T between and 5.5 T, includ- 
ing the critical field H c i — 2.02 T. Arrows mark the phase 
transition into the AFM ordered state. The inset displays 
in addition to a c , the coefficient a a {H || c, AL a ) and the 
calculated volumetric coefficient ay = 2a a + a c for T. 



C. Thermal Expansion 

The linear thermal expansion coefficient 

1 dALi 



L dT 



(13) 



is defined as the temperature derivative of the length 
change ALi along a certain crystallographic direction i. 
In tetragonal systems, such as DTN, the volumetric ex- 
pansion can be calculated from the linear coefficients 



ay = 2a a 



(14) 



along the crystallographic a and c direction. 

The main panel of Fig. [5] shows a c for H \\ c between 
and 5.5 T. For measurements well above H c i, the tran- 
sition into the ordered phase is indicated by a distinct 
anomaly, marked with arrows. In the temperature range 
up to 5 K no significant contribution from the lattice is 
observed. In zero field, a c shows a rather broad minimum 
that shifts to lower temperatures in higher fields and can 
be attributed to the thermal population of energetically 
higher spin states \S Z >— ±1, similar to the maximum 
in the specific heat. In zero field these spin states are 
equally occupied, because they have same energy. The 
largest AFM exchange along the c-axis leads to the dom- 
inant magnetostrictive effect. Since the PM ground state 
is a product of S z — state to a good approximation, 
the thermal excitation of S z = ±1 states increases the 
nearest-neighbor XY AFM correlations along the c-axis 
(S r • S r+ec ). This increase leads to an attractive mag- 
netostrictive force between nearest-neighbor ions along 



Thermodynamic quantity 


H cX (T) 


Magnetization M(H, T) 


2.08 


Magnetocaloric effect (MCE) 


2.06 


Magnetostriction AL || c 


2.02 


Magnetostriction AL _L c 


2.08 



TABLE II: Critical field H c i estimated by different experi- 
mental methods and setups. 



the c-axis that shrinks the lattice as the temperature in- 
creases from zero. On the other hand, the magnetrostric- 
tive force disappears at high enough temperature because 
(S r • S r +e c ) — > for T — > oo, implying that AL C must 
have a minimum at a finite temperature where a c changes 
sign. This expected behavior is fully consistent with the 
experimental results shown in Fig.[5j The critical field 
for the thermal expansion measurements was estimated 
by detailed magnetostriction measurements, whereas the 
magnetostriction coefficient 



A, 



1 dALj 
LL L dH 



(15) 



is definded as the magnetic field derivative of the length 
change along the i direction. The analysis of the data 
(not shown) gives H cl = 2.02 T for AL C and H cl = 2.08 T 
for AL a . The difference in the critical field values can be 
attributed to the application of small pressure on the 



sample during the experiment (see also section IV D ) . 

The inset of Fig. [5] compares a c , a a and ay for the 
zero-field measurement. Between 5K and 3K, the ther- 
mal expansion is dominated by the length change along 
the c direction, a c > 0, a a — 0, because of the quasi-1 
dimensional nature of the magnetic interactions in this 
temperature range. In the temperature range below 3 K, 
the thermal expansion coefficient a c is negative and has 
the minimum that is expected because the curve a c (T) 
must change sign at a finite temperature. In contrast, 
the thermal expansion coefficient in the plane a a is pos- 
itive for T < 3 K with a maximum at the temperature 
T = T m where the minimum occurs in a c . Calculating 
ay via Equation ( 14 ) reveals that the thermal expansion 
coefficients a a and a c strongly compensate each other 
and that the volume coefficient changes much less than 
the individual a, values. 



IV. DISCUSSION 

A. Critical Field H cl 

The values for the lower critical field H cl , estimated 
individually for the different experimental quantities, are 
summarized in Table [TTJ. The values vary by ± 30 mT 
around 2.05 T. This difference can be attributed to the 
application of small pressure on the sample during the ex- 
periment. The spring- loaded capacitive dilatometer for 
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FIG. 6: (color online) Experimental magnetization AM(T) 
versus temperature T (symbols) at the critical field hqH = 
2.08 T applied along the crystallographic c direction of DTN. 
The dashed line represents QMC results, whereas the solid 
line indicates the analytic calculations. The symbols in the 
inset show the magnetic Griineisen parameter T mag estimated 
from the data of the main panel and the specific heat values 
shown in Fig. [4] Dashed and solid lines represent QMC and 
analytical results respectively. 



thermal expansion and magnetostriction measurements 
can intrinsically apply a small amount of uniaxial pres- 
sure along the measured sample length. The misalign- 
ment of the crystal orientation out of H || c indicated 
by the magnetization experiments, is another reason for 
variations in H c \. 



B. Comparison Theory - Experiment 

Fig. [6] shows a comparison between the experimental 
magnetization (symbols) at the critical field and analyti- 
cal calculations (solid line) and QMC simulations (dashed 
line). We observe that the QMC and analytic results 
agree with the experimental data within the error bars 
below 0.5 K. A T 3 / 2 behavior is expected below 0.3 K for 
the BEC universality class, in contrast to T 2 dependence 
of an Ising-like QCP. While the analytic calculation is 
only valid at low temperatures (or low density of bosons) 
the QMC results remain valid at any temperature. This 
is the reason why the QMC results are in very good agree- 
ment with the experimental data up to 2.2 K. 

Figure[7] shows the experimental data of the specific 
heat, C mag /T(T), (symbols) at and slightly below the 
critical field H cl , compared with data of analytic calcu- 
lations (broken lines) and QMC simulations (solid line) 
in a double logarithmic display. The experimental data 
exhibit a VT dependence in the low temperature limit 
at the critical field - in agreement with mean-field cal- 
culations and QMC simulations - which is characteristic 



FIG. 7: (color online) Magnetic specific heat as C m a 3 /T as a 
function of temperature T at the critical field H c \ = 2.06 T 
(stars), at 2.04 T (open circles), 2 T (filled triangles) and 1.9 T 
(open diamonds) compared with QMC results at H c i (solid 
line) and with analytical calculations (dashed, dotted lines) 
at H c \ and for fields below the critical field. 



for the 3-dimensional BEC universality class. We observe 
that the experimental data agree well within error bars 
with the analytic results down to the lowest tempera- 
tures, whereas there is a slight deviation of the QMC 
results below 0.3 K. This discrepancy is most likely due 
to a small error in the numerical determination of the 
critical field H c ±. A linear temperature dependence of 
the specific heat C mag /T(T) is expected for an Ising-like 
QCP, which we can exclude from our measurements. The 
broad Schottky maximum around 2.5 K in the experimen- 
tal data originating from the population of \S Z >= ± 1 
excited spin states, can only be reproduced by the QMC 
simulations because they remain valid up to arbitrarily 
high temperatures. Seen in the experimental data as well 
as in the QMC results is the crossover from 3 dimensional 
(C/T ~ T 1 ' 2 ) to 1 dimensional (C/T ~ T" 1 / 2 ) behavior 
of the specific heat at higher temperatures for the mea- 
surement at H c \. The change of slope is marked with a 
broad hump around 0.5 K. Note that the 1-dimensional 
temperature dependence overlaps with the above men- 
tioned Schottky contribution to the specific heat. 

For H < H c i measurements, the analytical curves as 
well as the experimental values lie below the specific heat 
data at the critical field. They grow smaller as the dis- 
tance from H c \ increases. This observation confirms the 
correct estimation of the value H c \ = 2.06 T for the 
critical field in the MCE measurements. Furthermore, 
Fig. [7] shows nice agreement between analytical predic- 
tions and experimental data at low temperatures for all 
shown fields. Deviations seem to be larger for H < H c ±, 
but this is an effect of the double-logarithmic display. 

In presence of Ising-like anisotropy, the gap should re- 
open inside the AFM phase. We do not observe any ex- 
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FIG. 8: (color online) The main panel shows the normalized 
length change AL/L (symbols) as a function of the tempera- 
ture T of DTN for H || c ]| L c at the critical field H cl = 2.02 T 
in a semi-logarithmic display. The experimental data are 
compared with QMC calculations, where the spin-spin cor- 
relator (SSC) (dashed line) and the SSC together with the 
(S z ) 2 term (dotted line) are taken into account. The in- 
set (a) shows the experimental data at zero field (solid), the 
critical field H c i = 2.02 T (dotted line) and the difference 
AL c /L(H cl ) - AL c /L(0)) (symbols). In inset (b) the abso- 
lute values AL/L (symbols) are plotted as a function of T in 
comparison with ~ T 1/2 (solid), ~ T 3/2 (dashed), and ~ T 2 
(dotted line) temperature dependence. 



ponential temperature dependence in the measurements 
for fields above H c \ in Fig. gj namely 2.09 T, 2.14 and 
2.2 T. This, however, could also be due to (i) lack of data 
at temperatures below 80 mK and (ii) the onset of the 
phase transition seen as a broad anomaly in the specific 
heat. 

Before we discuss the scaling behavior of the thermal 
expansion coefficient at the critical field, let us have a 
closer look at the length change AL C /L, which is shown 
in Fig. [8] along with results from QMC simulations. The 
inset (a) of Fig. [8] shows the experimental length at the 
critical field after the subtraction of the data in 0T, in 
order to separate quantum critical from the non-critical 
magnetic contributions of the sample. The same pro- 
cedure was done for QMC data. In QMC simulations, 
AL C /L at H c i is obtained from the estimation of the spin- 
spin correlator (SSC), (S r • S r+eiy and optionally ad- 
ditional terms. The main panel of Fig. [8] shows the com- 
parison between the experimental data and QMC simu- 
lations. The qualitative features are well reproduced by 
the expectation value of the SSC for temperatures above 
0.3 K. The scaling factor between experiment and theory 
is 1.85 • 10~ 4 . AL C /L follows a T 1 power law with 7 be- 
tween 2 and 3/2 as illustrated in inset (b) of Fig. [8] This 
result is in close accordance with the expected BEC be- 
havior of ~ T 3 / 2 The discrepancy between experimental 



FIG. 9: The volume thermal expansion coefficient divided by 
T is presented for fields close to the critical fields as a function 
of temperature T in a semi-logarithmic plot for DTN with 
field direction H | c. The inset shows the absolute value of 
the thermal Gruneisen parameter T t h = ay /C at the critical 
field Hd- The solid line in the inset illustrates T -1 behavior. 



data and QMC simulations below 0.3 K can be attributed 
to additional contributions to AL C /L besides the SSC. 
However, an additional consideration of a ((S*) 2 } term 
that is expected from symmetry arguments, does not im- 
prove the agreement significantly. Therefore, the origin 
of this discrepancy remains unclear at the present level 
of analysis. 

Figure[9] summarizes the volumetric thermal expansion 
coefficient divided by the temperature ay /T for fields 
at and close to the QCP. Note, that according to equa- 
tion ( 14 ) different critical values H c i were considered for 
a„ and a r and taken into account for the estimation of 



ay. We find, that at H c \ the values ay /T have a weak, 
but finite temperature dependence (solid line) meaning 
that the thermal expansion coefficient ay diverges as ex- 
pected at the QCP. For magnetic fields H < H c i the 
low-temperature values ay /T show similar behavior. 

In general, pressure p and magnetic field H are equiva- 
lent parameters of the free energy F(p, H) in close vicin- 
ity to the critical field \H c i — H\ <C H c \. Therefore, ther- 
modynamic quantities derived from pressure and field de- 
pendencies can be converted 



dp dH 



(16) 



by multiplication with the prefactor 57 = dH c \ / dp, which 
is the hydrostatic pressure dependence of the critical field 
H c i. 



Equation ( 16 1 implies that the compressibility k = 
d 2 F/dp 2 is proportional to the magnetic susceptibility 
that is a step function of the magnetic field at H c i 
and is shown in Fig. 10 Therefore, k increases rapidly 
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FIG. 10: Volume magnetostriction coefficient Ay and suscep- 
tibility x = H^BM/dH at 0.1K of DTN. The susceptibility 
is scaled by the hydrostatic pressure dependence of the critical 
field Q. estimated from the Ehrenfest relations. 



at H c i and leads to a softening of the crystal lattice, 
recently demonstrated experimentally on DTN by ul- 
trasound experiments.^ Considering this, we speculate, 
that the huge change in the lattice properties is respon- 
sible for the deviation of ay from the exact behavior of 
a BEC. Similar crystal softening as a precursor for quan- 
tum criticality was e.g. also observed at the metamag- 
netic transition in CeRu2Si2p^ 

The thermal expansion coefficient ay for fields above 
H c i shows the clear onset of the phase transition with 
positive values cvy > 0. 



C. Griineisen Parameter 

The magnetic Griineisen parameter, T mag , is given in 
the inset of Fig. [6] and compared with theoretical data 
of the QMC simulation and the analytical calculations. 
In the temperature range below 0.3 K, the experimental 



mag 



shows the onset of divergence, as expected for a 
QCP. As far as it is observable in the low temperature 
limit, T mag follows the theoretical prediction ~ T -1 , be- 
cause the magnetization and specific heat obey the ex- 
pected behavior individually. For the thermal Griineisen 
parameter, r t ^, shown in the inset of Fig.[9j a T -1 behav- 
ior is found as well for temperatures below 0.6 K. These 
observations are in fully agreement with the universality 
class of a BEC QCP. 



M, allow us to obtain the pressure dependence of (i) the 
transition temperature T/v 



dT 



N 



dpi 



V m T, 



N 



AC n 



and of (ii) the critical field H cl 
dH cl 



dpi 



V r , 



(17) 



(18) 



by the Ehrenfest relations, which hold at phase tran- 
sitions of second order. Equations ( 17|18 ) are valid for 
uniaxial as well as hydrostatic pressure, dependent if the 
linear or volume coefficients A^ or a* are used in the anal- 
ysis. 

We find dT N /dp c = 18.5 K GPa -1 at H = 2.2 T for 
uniaxial pressure applied along the crystallographic c di- 
rection. This value is rather high, even compared to 
TICUCI3, which already exhibits a huge value of sev- 
eral Kelvin per GPa dependent on the considered field 
ranged Direct measurements of the thermal expansion 
AL a are currently not available, however, from measure- 
ments in small magnetic fields H < H c \ , we expect a neg- 
ative uniaxial pressure dependence 9Tjy /dp a that should 
partially cancel the value (18.5 K GPa -1 ) along c, lead- 
ing to a smaller but positive hydrostatic value. Similar 
behavior was observed in TlCuCFj!^ 

The pressure dependence of the critical field, 
dH c i/dpi, could be obtained at 0.1K for both directions 
c and a. Because of a negative jump AA C in the mag- 
netostriction coefficient at H c i (data not shown), the 
value dH c i/ dp c = - 6.76 T GPa -1 is negative. In con- 
trast, 9-ff c i/<9p a =1.84TGPa -1 is positive and the result- 
ing hydrostatic pressure depe nde nce, which is equivalent 
to the prefactor in equation 16 adds to - 3.07T GPa -1 . 
The correct estimate of f2 can be proved by a compar- 
ison of Ay with the magnetic susceptibility \ measured 
at 0.1K. Both thermodynamic quantities are linked via 



equation ( 16 ) and it follows 



n 



d 2 F 



dH 2 



(19) 



d 2 F 
dpdH 

FigureJTU] shows an excellent agreement between the mag- 
netostriction coefficient Ay and the scaled susceptibility 
ilX- We conclude from the analysis of the Ehrenfest re- 
lations, that the application of uniaxial pressure p c along 
the c direction increases the ordering temperature Tjv 
and reduces the critical field H&. The behavior is op- 
posite for uniaxial pressure along the a direction. The 
response of DTN to hydrostatic pressure is dominated by 
the uniaxial pressure dependence along the c axis because 
this is the direction of the dominant magnetic exchange 
interaction. 



D. Ehrenfest relations 

The anomalies at the phase boundaries observed in the 
thermal expansion coefficient, c^, the specific heat, C, the 
magnetostriction coefficient, Ai, and the magnetization, 



V. SUMMARY 

We present a comprehensive experimental and theo- 
retical study of the thermodynamic properties: specific 
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heat, magnetization, and thermal expansion in the 
vicinity of the field-induced QCP around H c \ w 2 T 
in NiCl2-4SC(NH 2 )2- This point marks the entrance 
into an 3-dimensional XY antifcrromagnctically ordered 
state and can be described within the formalism of a 
BEC of magnons. We find a T 3 / 2 low temperature 
behavior of the specific heat and the magnetization at 
H c i that are in agreement with the universality class of 
a BEC-QCP. The thermal expansion coefficient shows 
a temperature dependence T 7 with 3/2 < 7 < 2 for 
T — > 0, which is in close agreement with expectations 
for this kind of QCP. QMC simulations nicely reproduce 
the features observed in the magnetization, specific heat 
and thermal expansion over a broad temperature range. 
Only the low temperature dependence in the specific 
heat and thermal expansion deviates due to intrinsic 
uncertainties in the simulation method. Furthermore, 
we analysed the thermal, T t h, as well as the magnetic 
Griineisen parameter T magi which are key quantities 
for the identification of QCPs and diverge with specific 
power laws. Experimentally, we found a T" 1 diver- 
gence for T mag and T th as expected for a BEC-QCP. 
Moreover, we estimated the influence of pressure on the 
transition temperature, Tjv, and the critical field, H c \ 



and found opposite effects for uniaxial pressure along 
the crystallographic a and c axes. Uniaxial pressure 
along c (a) direction increases (reduces) the ordering 
temperature Tjv while it reduces (increases) the critical 
field H c \. Due to the quasi- 1 dimensional character of 
the exchange interactions in DTN, the application of 
hydrostatic pressure increases T N and reduces H cl . Our 
results encourage pressure experiments that shift H cl 
to zero field. Since the field induced QCP is BEC-likc, 
the pressure induced QCP should belong to the 0(2) 
universality class in dimension D = 3 + 2. 
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